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OPTIMAL,  VIBRATION  REDUCTION  OVER  A  FREQUENCY  RANGE 


Abstract*  This  is  a  review  of  optimal  vibration  reduction 
techniques  for  systems  subject  to  harmonic  excitation  over  a 
frequency  range.  Only  passive  means  of  control  are  considered.  The 
objective  functions  used  for  optimization  are  restricted  to  those 
which  relate  directly  to  some  measure  of  frequency  response.  Other 
common  optimization  goals  such  as  weight  minimization  with  response 
constraints  are  not  included  in  this  survey. 


CONVENTIONAL  DYNAMIC  VIBRATION  ABSORBER 


The  simplest  device  used  to  attenuate  the  steady-state 
vibration  of  a  mechanical  system  over  a  frequency  range  is  the 
conventional  dynamic  absorber.  In  the  classical  analysis , the  main 
system  is  modelled  as  a  mass  resonating  on  a  spring  (Figure  1).  The 
absorber  is  tuned  to  resonate  such  that  the  motion  if  its  mass  m 
becomes  relatively  large  and  the  motion  of  the  main  mass  m  is  2 
minimized.  The  first  analysis  of  the  absorber  is  usually  attributed 
to  Ormondroyd  and  Den  Hartog  [1],  A  detailed  discussion  of  optimal 
tuning  and  damping  parameters  is  given  in  Den  Hartog' s  book  [2].  To 
write  the  equations  of  motion  in  a  dimensionless  form,  the 
following  symbols  are  introduced 
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static  displacement  of  main  mass  produced  by  a 
force  F 

mass  ratio  «  absorber  mass /main  mass 
uncoupled  natural  frequency  of  main  system 
uncoupled  natural  frequency  of  damper  system 
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ratio  of  natural  frequencies 

ratio  of  the  exciting  frequency  to  the 

uncoupled  natural  frequency  of  main  system 

critical  damping 

damping  ratio 

1,2)  -  displacement  ratio 


With  this  notation  the  equations  of  motion  in  the  frequency 
domain  |  may  be  written 
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If  A/C  »  B/D,  then  z  becomes  equal  to  A/C  which  is  independent  of 

damping  £  .  Por  fixed  m  and  f ,  the  condition  A/C  -  B/D  yields  a 

2 

quadratic  equation  in  g  the  two  roots  of  which  are  the  frequencies 
^l'g2  w^ere  t*ie  amplitude  I  I  of  the  main  mass  is  independent  of 
damping.  The  points  ( g^ , I 9^ ) 1  )  and  ( g2 ,  I  z1<g2 )  l  )  on  the  frequency 
response  curve  lz1(g)l  are  called  fixed  points  or  invariant  points. 
Den  Hartog  optimizes  absorber  performance  by  first  choosing  f  -  n 
(tuning)  so  that  the  fixed  points  are  adjusted  to  equal  height  and 
then  finding  the  value  of  £  (optimum  damping)  to  make  the  frequency 
response  curve  pass  through  one  of  the  fixed  points  with  a  horizontal 
tangent.  His  result  for  optimum  tuning  is 


To  find  optimum  damping  for  fixed  substitute  equation  (5)  into 
equation  (3),  find  the  derivative  of  I  z^  I  with  respect  to  gf  then 
evaluate  at  one  of  the  fixed  points,  say  g^ ,  and  equate  the  result 
to  zero.  From  the  equation  thus  obtained,  the  damping  ratio  may  be 
found  a9 
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On  the  other  hand,  if  the  fixed  point  chosen  is  the  one  at  a^,  the 
resulting  damping  ratio  is  given  by 
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Den  Hartog  recommends  the  use  of  the  average  value 
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EXTENSIONS  OF  THE  CONVENTIONAL  DYNAMIC  ABSORBER 


Lewis  [3]  extended  this  optimization  procedure  to  multiple 
dearee-of-freedom,  discrete,  undamped  systems  to  which  a 
conventional  viscous  vibration  absorber  is  attached.  For  such  a 
system  with  N  masses,  Lewis  proves  the  existence  of  2N-2  fixed 
points.  Thus,  if  the  designer  selects  a  resonant  peak  of  primary 
interest,  the  adjacent  invariant  points  in  the  vicinity  of  this  peak 
can  be  adjusted  as  described  above  to  obtain  optimum  tuning  and 
damping.  As  a  particular  case,  Lewis  analyzes  the  two-mass  system 
shown  in  Figure  2  and  shows  that  the  invariant  points  occur  at  the 
frequencies 
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with  optimum  tuning  given  by 
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and  optimum  damping  by 
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Snowdon  [4]  modified  the  conventional  absorber  by  adding  a  spring 
in  series  with  the  damper  and  showed  that  the  resulting  device,  called 
a  three-element  dynamic  absorber,  active  through  a  greater  frequency 
range,  reduced  transmissibility  (absolute  value  of  the  ratio  of  main 
mass  displacement  to  the  prescribed  displacement  of  a  foundation 
attached  to  the  main  mass  by  a  spring  )  more  effectively  at  the  center 
of  this  range.  In  the  same  paper,  Snowdon  introduced  dual  absorbers 
to  create  a  pronounced  trough  in  the  transmissibility  curve  while 
avoiding  the  two  large  compensating  resonant  peaks  that  appear  in  the 
transmissibility  curve  of  the  single  absorber.  Snowdon  attaches  these 
absorbers  to  an  undamped  single  degree-of- freedom  system  and  provides 
design  information  for  them  in  graphical  form.  Two  of  his  curves  have 
been  reproduced  here  to  show  the  changes  in  transmissibility  (Figures 
3  and  4). 


Randall  et.  al.  [5]  provided  optimum  absorber  design  curves  for 
primary  systems  with  deunping,  (Figure  5 ).  They  proved  that  once 
damping  is  introduced  into  the  main  system,  invariant  points  no 
longer  exist.  Defining  a  performance  measure  G  by 


G  -  max  I x, (w) I 
u>  1 


(12) 


they  used  a  numerical  search  to  compute  tuning  and  damping  parameters 
that  minimize  G.  Soom  and  Lee  [6]  solved  the  same  problem  by 
nonlinear  programming  and  investigated  the  possibility  of  using 
nonlinear  springs  to  improve  broad  band  response.  They  also  tried 
several  different  objective  functions,  for  example 

G.  =  E(lx  (cj)  -  11)^  for  u  such  that  lx  ( )  I  >  1 
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Nonlinearities  have  also  been  treated  by  Roberson  [7]  and  Arnold 
[8).  Roberson  considered  the  system  shown  in  Figure  6  where  the 
spring  between  the  masses  and  is  nonlinear.  Its  load- 
deflection  curve  is  the  sum  of  a  linear  and  a  cubic  term.  Since  the 
'^‘■'orber  is  undamped  in  this  case,  its  effectiveness  is  limited  to  a 
small  frequency  range.  Roberson's  synthesis  criterion  is  the 
maximization  of  this  suppression  band.  Arnold  [8]  studied  the  same 
system  as  Roberson  but  set  c  -  0.  He  provides  frquency  response 
curves  of  the  system  for  hardening  and  softening  coupling  springs  as 
well  as  for  a  softening  spring  designed  according  to  Roberson's 
optimum  system-parameter  specifications.  In  both  [7]  and  [8]  the 
restoring  force  in  the  nonlinear  spring  as  a  function  of  spring 
extension  is  taken  of  the  form 

R(x)  -  c(x  ±  m2x3)  (14) 

where  the  plus  sign  indicates  a  hardening  spring  characteristic  and 
the  minus  sign  indicates  a  softening  one. 

KwaX  et.  al.  [9]  and  Haug  and  Arora  [101  used  a  steepest 
descent  algorithm  to  solve  the  classical  absorber  problem  over  a 
finite  frequency  range.  Their  objective  function  is 
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where  z is  as  defined  in  equation  (1)  and  the  frequency  interval 
of  interest  is  [a,b].  They  minimize  G  subject  to  the  constraints 
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where  the  notation  is  as  in  equation  ( 1 ) ,  They  showed  that  for 
certain  finite  frequency  intervals,  designs  superior  to  Den  Hartog ’s 
infinite  frequency  range  optimum  may  be  found. 

SYSTEMS  WITH  CONTINUOUS  MEMBERS 


Plunkett  [11]  proved  the  existence  of  invariant  points  for 
undamped  continuous  systems  to  which  a  single  discrete  damper 
has  been  connected  and  exploited  this  property  to  determine  optimum 
damping  for  these  systems.  His  approach  is  a  generalization  of  the 
methods  of  Den  Hartog  [1,2]  and  Lewis  [3].  He  considers  the 
vibration  velocity  at  one  point  of  a  linear  system  resulting  from  a 
sinusoidal  force  at  another  point.  Let  the  force  applied  at  a  point 
1  be  Fj  and  v^  be  the  velocity  at  the  same  point.  Suppose  a  damper 
with  damping  constant  c  is  applied  across  points  with  relative 
velocity  v^  and  let  the  unknown  vibration  velocity  at  the  point  of 
interest  be  .  Plunkett  shows  that 
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and  notes  that  equation  (19)  has  the  same  form  as  Den  Hartog’ s 
equation  (3).  Thus  when  *>2/^3  "  the  ratio  y$/Fi  is  independent 
of  c  and  has  the  invariant  points  property  described  above  in 
conjunction  with  the  conventional  dynamic  vibration  absorber.  The 
value  of  c  that  gives  a  zero  slope  of  the  amplitude  Iv^F^I  with 
respect  to  frequency  (J  at  the  invariant  value  of  Iv^/F^l  is  an 
optimum  value  for  c.  Henney  and  Raney  [12]  applied  a  similar 
technique  to  optimize  damping  for  vibrating  uniform  beams  and 
studied  the  sensitivity  of  maximum  displacement  response  to 
deviation  from  optimum  damping. 

Snowdon  [13]  ,[14]  described  how  absorbers  could  reduce  the 
force  transmitted  to  the  terminations  of  undamped  cantilever  beams  at 
the  resonant  frequencies  of  beam  vibration.  The  systems  under 
consideration  are  shown  in  Figure  7.  The  force  transmissibility  T 
across  the  beam  is  defined  as  follows 
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where  F^  is  the  force  transmitted  to  the  fixed  end  and  is  the  applied 
force.  When  a  single  absorber  is  used,  as  in  Figures  7a  and  7b. 
essentially  the  same  procedure  as  that  used  for  the  conventional 
absorber  is  applicable.  The  fixed  points  are  located  by  comparing 
transmissibility  curves  obtained  for  zero  and  infinite  damping. 

The  frequency  ratio  is  then  chosen  so  that  the  fixed  points  lie  on 
transmissibility  curve  actually  take  the  equal  values  of 
transmissibility  at  the  fixed  points.  When  two  absorbers  are  applied 
simultaneously  to  the  beam.  Figure  7c.  Snowdon’s  approach  i3  to 
assign  each  absorber  the  values  of  optimum  tuning  and  damping  that 
were  determined  when  the  absorber  Weis  attached,  individually,  to  its 
present  position  on  the  beam.  Optimum  supression  of  both  the  first 
and  the  second  beam  resonances  results  from  this  procedure. 

Jacquot  and  Foster  [ 15 ]  considered  an  undamped  single-degree- 
of- freedom  system  equipped  with  a  damped  cantilever  beam  serving  as 
a  vibration  absorber.  The  authors  developed  an  approximate  system 
dynamics  by  an  assumed  modes  approach  retaining  only  one  mode  of  the 
beam.  Since  the  resulting  equations  were  of  the  form  obtained  for 
the  conventional  absorber  problem,  the  same  tuning  and  damping 
approach  was  applicable.  Snowdon  [16]  analyzed  platelike  absorbers 
attached  either  to  a  lumped  mass-spring  system  or  to  a  plate  with 
small  internal  damping  and  presented  optimum  design  parameters  in 
graphical  form. 

Jacquot  [17]  appended  a  discrete  absorber  with  structural  or 
viscous  damping  to  a  beam.  Figure  8.  Using  an  assumed  mode  approach 
as  in  [15],  he  derived  an  approximate  frequency  response  given  by 

i/rfVA)  =*  T2  -  A2 _  _  (21) 
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and  0.  ,  f>  are  the  beam  mode  shapes  and  eigenvalues  satisfying 
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Because  the  form  of  equation  ( 21 >  i9  the  same  as  that  found  in  the 
conventional  absorber  problem  with  the  exception  that  /x  is  replaced 
by  tip  (a>,  Jacquot  was  able  to  tune  and  damp  to  flatten  the 
fundamental  resonance  response  using  Den  Hartog's  procedure. 

Warburton  and  Ayorinde  [18]  extended  Jacquot’s  representative 
mode  method  to  plates  and  shells  by  defining  appropriate  effective 
mass  and  stiffness  for  such  elastic  bodies.  The  accuracy  of  this 
single-mode  approximation  was  favorably  affected  if  adjacent 
resonant  frequencies  were  well  separated  from  the  natural  frequency 
for  which  the  absorber  was  being  tuned  and  adversely  affected  as  the 
absorber  size  increased.  In  a  companion  paper  fl9),  the  authors 
considered  cylindrical  shells  as  examples  of  dynamically  complex 
structures,  for  which  the  ratio  of  adjacent  natural  frequencies 
tends  toward  unity.  They  showed  that  as  dynamic  complexity 
increased  optimum  absorber  parameters  deviated  from  those  calculated 
for  an  equivalent  single-dearee-of-freedom  system. 


MULT  I  PARAMETER  DESIGN  BY  NONLINEAR  PROGRAMMING 


McMunn  and  Plunkett  [20]  developed  a  computational  method  to 
optimize  multiple  dampers  for  large  mechanical  systems.  Damping  was 
defined  to  be  optimum  if  the  maximum  response  over  a  range  of 
excitation  frequencies  was  minimized.  The  response  function  of 
interest .  for  example 


f(w,c)  - 


where  x.  is  the  displacement  of  the  ith  mass  and  P.  is  the  force  on 
the  jth  mass,  is  first  maximized  with  respect  to  excitation 
frequency  (J  using 


and  then  minimized  with  respect  to  the  vector  c  of  damping  values 


f (w, c) 


In  fact,  the  minimization  over  c  also  considers  values  of  f 
at  w  •  0  ,  and  w  ■  ft  if  a  finite  frequency  interval  is  *>*lng 

considered.  Two  multlple-degree-of -freedom,  multiple-damper 
discrete  systems  and  a  column  with  complex  modulus  damping  were 
studied  as  examples  using  this  approach.  A  similar  approach  was 
adopted  by  Ng  and  Cunniff  [21]  who  designed  a  three 
degree-of-freedom  isolation  system  and  verified  their  results  by 


experimental  tests.  In  their  formulation,  the  authors  define  a 
primitive  function 


2  -1  -1  * 
1^1  -  (P  f)(P  f) 


(27) 


where  *  denotes  complex  conjugate  transpose, 
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and  f  is  the  force  vector.  Their  first  step  is  to  minimize 

over  damping  variables  the  primitive  function  maximized  over  frequency 
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They  repeat  this  procedure  until  an  optimum  is  reached 
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Dale  and  Cohen  [22]  extended  the  method  of  McMunn  and  plunkett 
to  continuous  systems  whose  steady-state  equations  of  motion  could  be 
reduced  to  a  set  of  ordinary  differential  equations  containing 
spatial  coordinates  as  independent  variables.  Both  dissipative  and 
nondissipative  design  parameters  were  included. 


Lunden  [23]  presented  a  nonlinear  programming  solution  based  on 
a  sequential  unconstrained  minimization  technique  to  the  problem  of 
determining  a  continuous  damping  distribution  which  minimizes  the 
maximum  response  of  a  vibrating  beam  over  a  specified  frequency 
interval.  In  a  second  paper  [24],  the  author  applied  the  same 
approach  to  vibrating  frames.  In  these  references,  the  maximum 
response  F  in  the  frequency  interval  studied  is  written  as 
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where  T)  denotes  distributed  structural  damping  and  i?8  denotes  the 
structural  damping  constants  for  discrete  springs  in  the 
system.  An  exact  displacement  method  is  used  with  hysteretic 
damping  introduced  by  the  loss  factor  tj  giving  a  complex  bending 
stiffness  EI(  1  +  The  4x4  stiffness  matrix  for  a  beam 

element  then  takes  the  form 


E  -  El  ( 1  +  jT)J)K  x 


(32) 


where  K  is  a  4  x  4  matrix  of  transcendental  functions  of  frequency. 


Kitis  [25]  utilizes  structural  reanalysi9  and  modal  techniques 
with  nonlinear  programming  to  make  tractable  problems  in  which  the 
systems  under  consideration  contain  a  large  number  of  degrees  of 
freedom.  The  repetitive  computations  of  response  required  in  the 
nonlinear  programming  portion  of  the  optimal  design  are  carried  out 
using  efficient  reanalysis  methods  or  condensed  eigenproblem 
solutions  so  that  computation  time  in  the  structural  analysis  phase 
of  the  design  is  reduced. 

IMPEDANCE  MATCHING 

The  progressive  wave  solution  to  the  wave  equation  for 
continuous  chain-like  systems  has  been  used  to  minimize  vibratory 
response  over  a  frequency  range  by  means  of  "impedance  matching". 

To  illustrate  the  basic  idea  of  this  technique  consider  the  rod 
shown  in  Figure  9.  [14].  Designating  the  displacement  at  any  point 
x  by  £ ,  the  wave  equation  may  be  written  as 
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If  the  end  of  the  rod  at  x  ■  o  is  subjected  to  a  sinusoidally 
varying  force  P  .  then  the  driving  point  impendance  is  given  by 


Z„  = 


_Fo _ 

onidt)xm0 


Fo 

jui£(x-o) 


(35) 


Snowdon  [14]  shows  that  the  impedance  Z  defined  in  equation 
( 35 ) ,  may  be  written  as 


/  -  ReJtf 
/  +  ReJ<f 


(36) 


where  Re^P  describes  the  relative  magnitude  of  and  the  phase 
difference  between  the  incident  and  the  reflected  waves.  The 
characteristic  impedance  Z  is  defined  as  the  impedance  of  an 
infinitely  long  rod  in  which  reflections  do  not  occur,  the  value  of 

Z  ±b  obtained  from  eauation  ( 36 )  by  equatina  R  to  zero 
ch 


Z 


ch 


(37) 


A  matched  condition  will  occur  if  a  damper  of  dampina  constant  Z 

*  *  ch 

is  attached  to  the  rod  at  x  -  0.  Thi9  attachment  causes  the  ratio  R 
of  reflected  to  incident  wave  to  be  zero  so  that  no  vibration 
response  buildup  due  to  reflected  waves  is  possible.  This  idea  has 
been  applied  to  shafts  on  supports  [26].  Here  the  dynamic  response  is 
expressed  in  progressive  wave  form  like  electrical  response  waves  in 
tran9misison  line  theory.  Waveforms  for  a  uniform  shaft  flexibly 
supported  on  two  rotational  and  translational  mass-spring-damper 
units  at  the  ends  and  one  such  unit  in  the  interior  are  obtained. 

The  terminating  impedance  is  made  equal  to  the  characteristic 
impedance  of  the  shaft  to  obtain  a  matched  condition. 


LIMITING  PERFORMANCE  METHODS 


A  limiting  performance  approach  was  applied  in  Refs.  [27]  and 
[28]  to  the  optimal  design  of  vibratory  systems  over  a  frequency 
range.  In  this  approach,  instead  of  fixing  the  design  configuration 
at  the  outset,  those  parts  of  the  system  to  be  designed  are  replaced 
by  control  forces.  Then,  for  a  selected  cost  function  and  design 
constraint  the  absolute  optimal  performance  of  the  system  is 
computed  by  solving  an  optimization  problem  in  which  the  control 
forces  are  unknowns.  The  solution  is  called  the  limiting 
performance  of  the  system.  After  the  limiting  performance 
characteristics  have  been  found,  the  designer  can  choose  a 
prospective  configuration  for  the  part  of  the  system  to  be  designed 
and  apply  parameter  identification  techniques  to  find  optimum  design 
variable  values  so  that  the  designed  system  responds  as  closely  as 
possible  to  the  limiting  performance  response.  The  limiting 
performance  characteristics  are  found  by  linear  programming  and 
parameter  identification  can  be  accomplished  by  such  curve  fitting 
techniques  as  least  squares.  This  two-stage  procedure  has  been 
demonstrated  for  the  optimal  design  of  rotor  suspension  systems  in 
[27]. 


BOOKS  AND  MONOGRAPHS 


The  book  by  Haug  and  Arora  [10]  contains  considerable  material 
applicable  to  frequency  response  shaping,  although  specialized 
aspects  are  not  treated  in  detail.  Two  other  useful  references  are 
the  monograph  by  Sevin  and  Pilkey  [28]  and  the  boc>  by  Snowdon  [14]. 
Sevin  and  Pilkey  present  an  introduction  to  the  subject  and  a 
summary  of  the  state-of-the  art  up  to  1971.  Snowdon* s  book  contains 
a  wealth  of  information  on  vibration  absorbers  and  reduction  of  beam 
vibrations . 
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system,  =  mass  of  main  system,  =  mass  of 
absorber,  w  =  excitation  frequency. 

Figure  3.  Transmissibility  of  the  Three-Element 
Dynamic  Absorber 
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Figure  4.  Transmissibility  of  the  Dual 
Absorber 


Figure  5.  Vibration  Absorber  for 
a  Damped  System 


a.  Attached  at  end 


b.  Attached  at  midpoint 


c.  Attached  at  the  end  and  midpoint 

K  »  absorber  spring  constant 
a 

n  «  absorber  damping  constant 
a 

M  a  absorber  mass 
a 

Z  =  absorber  impedance 

Figure  7.  Dynamic  Absorbers  Attached  to  a  Cantilever 

Beam  Excited  by  a  Sinusoidally  Varying  Force 
at  its  Free  End 
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OPTIMAL  VIBRATION  REDUCTION  OVER  A  FREQUENCY  RANGE 


Abstract.  This  is  a  review  of  optimal  vibration  reduction 
techniques  for  systems  subject  to  harmonic  excitation  over  a 
frequency  range.  Only  passive  means  of  control  axe  considered.  The 
objective  functions  used  for  optimization  are  restricted  to  those 
which  relate  directly  to  some  measure  of  frequency  response.  Other 
common  optimization  goals  such  as  weight  minimization  with  response 
constraints  axe  not  included  in  this  survey. 

CONVENTIONAL  DYNAMIC  VIBRATION  ABSORBER 

The  simplest  device  used  to  attenuate  the  steady-state 
vibration  of  a  mechanical  system  over  a  frequency  range  is  the 
conventional  dynamic  absorber.  In  the  classical  analysis, the  main 
system  is  modelled  as  a  mass  resonating  on  a  spring  (Figure  1).  The 
absorber  is  tuned  to  resonate  such  that  the  motion  if  its  mass  m 

2 

becomes  relatively  large  and  the  motion  of  the  main  mass  is 
minimized.  The  first  analysis  of  the  absorber  is  usually  attributed 
to  Ormondroyd  and  Den  Hartog  [1].  A  detailed  discussion  of  optimal 
tuning  and  damping  parameters  is  given  in  Den  Hartog' s  book  [2].  To 
write  the  equations  of  motion  in  a  dimensionless  form,  the 
following  symbols  are  introduced 

x  ■  plk,  “  static  displacement  of  main  mass  produced  by  a 

S  v  1 

force  P 

m  m./m  *  mass  ratio  «  absorber  mass  /main  mass 

2  l 

ft  •  fk_An  )  -  uncouDled  natural  frequency  of  main  svstem 

n  11  2/2 

U)  «  (k  /ra  )  -  uncoupled  natural  frequency  of  damper  system 

n  2  2 

(1) 

f  »  w  /ft  ■  ratio  of  natural  frequencies 

n  n 

g  -  cj/ftn  "  ratio  exciting  frequency  to  the 

uncoupled  natural  frequency  of  main  system 

C  -  2m  n  -  critical  damDina 

c  2  n 

{  -  c2/Cc  "  damPin9  ratio 

z.  «  x.  lx  ( i  *  1,2)  m  disolacement  ratio 
1  i  st 

With  this  notation  the  equations  of  motion  in  the  frequency 
domain  may  be  written 

(1  -  g2)z1  +  /xf2(z1  -  z?)  +  ■)2nig(z1  -  z ^  -  1 

(2) 

-  g2z2  +  f2(z2  “  zi^  +  32g{(z2  -  z^)  -  0 


which  gives 


1  C  +  jeD 

where 


B  -  2g 


2  2  2  2  2 
C  -  -nf  g  +  (g4  -  l)(g4  -  f4) 

D  -  2g(l  -  g2  -  Ag2) 


(3) 


(4) 


If  A/C  -  B/D,  then  z becomes  equal  to  A/C  which  is  independent  of 
damping  For  fixed  /i  and  f ,  the  condition  A/C  -  B/D  yields  a 
quadratic  equation  in  g2  the  two  roots  of  which  are  the  frequencies 
gi'g2  where  the  amplitude  I  I  of  the  main  mass  is  independent  of 
damping.  The  points  (g^,  )  and  (g^,  )  on  the  frequency 

response  curve  \ z^( g)J  are  called  fixed  points  or  invariant  points. 
Den  Hartog  optimizes  absorber  performance  by  first  choosing  f  »  <j  / 
(tuning)  so  that  the  fixed  points  are  adjusted  to  equal  height  and 
then  finding  the  value  of  £  (optimum  damping)  to  make  the  frequency 
response  curve  pass  through  one  of  the  fixed  points  with  a  horizontal 
tangent.  His  result  for  optimum  tuning  is 


To  find  optimum  damping  for  fixed  M.  substitute  equation  (5)  into 
equation  (3),  find  the  derivative  of  I  z1 I  with  respect  to  g,  then 
evaluate  at  one  of  the  fixed  points,  say  g  ,  and  equate  the  result 
to  zero.  From  the  equation  thus  obtained,  the  damping  ratio  may  be 
found  as 


f 2  _  -  yluJji  +  2~> 

*  "  8(1  +  n)3 


(6) 


On  the  other  hand,  if  the  fixed  point  chosen  is  the  one  at  g7,  the 
resulting  damping  ratio  is  given  by 


ttf 3  ±  Mill U.  i  2n 


8(1  +  n)3 


(7) 


Den  Hartog  recommends  the  use  of  the  average  value 


l 


-ilL. 


8(1  +  m) 3 


(8) 


EXTENSIONS  OP  THE  CONVENTIONAL  DYNAMIC  ABSORBER 


Lewis  [3]  extended  this  optimization  procedure  to  multiple 
degree-of-freedom,  discrete ,  undamped  systems  to  which  a 
conventional  viscous  vibration  absorber  is  attached.  For  such  a 
system  with  N  masses*  Lewis  proves  the  existence  of  2N-2  fixed 
points.  Thus,  if  the  designer  selects  a  resonant  peak  of  primary 
interest,  the  adjacent  invariant  points  in  the  vicinity  of  this  peak 
can  be  adjusted  as  described  above  to  obtain  optimum  tuning  and 
damping.  As  a  particular  case,  Lewis  analyzes  the  two-mass  system 
shown  in  Figure  2  and  shows  that  the  invariant  points  occur  at  the 
frequencies 


2  K(M^  +  +  m)  KV2M1M2ro(M1  +  M?  +  m) 

‘Vb  "  M!(M2  +  m)  ±  2M1M2(M2  +  m) 


(9) 


with  optimum  tuning  given  by 


k 

k 


reT2M-ML  +  2M^  +  2M?m  -  M^m) 
2M1fM2  +  m)2 


(10) 


and  optimum  damping  by 


2 

c 


2  3 

K(8M ^  +  8Mp  +  SM^m  -  M^nQin 

4(M2  +  m)3(M1  +  M2  +  m) 


(11) 


Snowdon  [4]  modified  the  conventional  absorber  by  adding  a  spring 
in  series  with  the  damper  and  showed  that  the  resulting  device,  called 
a  three-element  dynamic  absorber,  active  through  a  greater  frequency 
range,  reduced  transraissibility  (absolute  value  of  the  ratio  of  main 
mass  displacement  to  the  prescribed  displacement  of  a  foundation 
attached  to  the  main  mass  by  a  spring  )  more  effectively  at  the  center 
of  this  range.  In  the  same  paper,  Snowdon  introduced  dual  absorbers 
to  create  a  pronounced  trough  in  the  transmissibility  curve  while 
avoiding  the  two  large  compensating  resonant  peaks  that  appear  in  the 
transmissibility  curve  of  the  single  absorber.  Snowdon  attaches  these 
absorbers  to  an  undamped  single  degree-of-freedom  system  and  provides 
design  information  for  them  in  graphical  form.  Two  of  his  curves  have 
been  reproduced  here  to  show  the  changes  in  transmissibility  (Figures 
3  and  4). 


Randall  et.  al.  [51  provided  optimum  absorber  design  curves  for 
primary  systems  with  damping,  (Figure  51.  They  proved  that  once 
damping  is  introduced  into  the  main  system,  invariant  points  no 
longer  exist.  Defining  a  performance  measure  G  by 


G  -  max  I  x.  (cj)  I 

oi  1 


(12) 


they  used  a  numerical  search  to  compute  tuning  and  damping  parameters 
that  minimize  G.  Soora  and  Lee  [6]  solved  the  same  problem  by 
nonlinear  programming  and  investigated  the  possibility  of  using 
nonlinear  springs  to  improve  broad  band  response.  They  also  tried 
several  different  objective  functions,  for  example 

G  -  E(Ix.(gj)  -  II)2  for  w  such  that  lx  (w)  I  >  1 
1  u>  1  1 


G_  -  max(wl x  (w) I ) 

2  0)  1 

G,  -  £lx. (w) I  2 

3  u  1 

G  -  £  cjI  X  (u)  I  2 

4  cj  1 


(13) 


Nonlinearities  have  also  been  treated  by  Roberson  [7]  and  Arnold 
[8).  Roberson  considered  the  system  shown  in  Figure  6  where  the 
spring  between  the  masses  and  is  nonlinear.  Its  load- 
deflection  curve  is  the  sum  of  a  linear  and  a  cubic  term.  Since  the 
absorber  is  undamped  in  this  case,  its  effectiveness  is  limited  to  a 
small  frequency  range.  Roberson’s  synthesis  criterion  is  the 
maximization  of  this  suppression  band.  Arnold  [8]  studied  the  same 
system  as  Roberson  but  set  c  -  0.  He  provides  frquency  response 
curves  of  the  system  for  hardening  and  softening  coupling  springs  as 
well  as  for  a  softening  spring  designed  according  to  Roberson's 
optimum  system-parameter  specifications.  In  both  [7]  and  [8]  the 
restoring  force  in  the  nonlinear  spring  as  a  function  of  spring 
extension  is  taken  of  the  form 

R(x )  -  c(x  ±  m2x3)  (14) 

where  the  plus  sign  indicates  a  hardening  spring  characteristic  and 
the  minus  sign  indicates  a  softening  one. 

Kwak  et.  al.  [9]  and  Haug  and  Arora  [10]  used  a  steepest 
descent  algorithm  to  solve  the  classical  absorber  problem  over  a 
finite  frequency  range.  Their  objective  function  is 

■  «?S*b  Ui(u)l 


G 


(15) 


where  z is  as  defined  in  equation  ( 1 )  and  the  frequency  interval 
of  interest  is  [a,b].  They  minimize  G  subject  to  the  constraints 


,?3*b 


-  X. 


f  <  f  C  f 
MIN  MAX 


;  <  p  <  p 

“MIN  *  “MAX 


where  the  notation  is  as  in  equation  ( 1 ) .  They  showed  that  for 
certain  finite  frequency  intervals,  designs  superior  to  Den  Hartog’s 
infinite  frequency  range  optimum  may  be  found. 

SYSTEMS  WITH  CONTINUOUS  MEMBERS 


Plunkett  [11]  proved  the  existence  of  invariant  points  for 
undamped  continuous  systems  to  which  a  single  discrete  damper 
has  been  connected  and  exploited  this  property  to  determine  optimum 
damping  for  these  systems.  His  approach  is  a  generalization  of  the 
methods  of  Den  Hartog  [1,2]  and  Lewis  [3],  He  considers  the 
vibration  velocity  at  one  point  of  a  linear  system  resulting  from  a 
sinusoidal  force  at  another  point.  Let  the  force  applied  at  a  point 
1  be  P1  and  be  the  velocity  at  the  same  point.  Suppose  a  damper 
with  damping  constant  c  is  applied  across  points  with  relative 
velocity  v  and  let  the  unknown  vibration  velocity  at  the  point  of 
interest  be  v^ .  Plunkett  shows  that 


*  =  ,b,  1  t-Jhii 

F,  J  l  +  jbjC 


and  notes  that  equation  (19)  has  the  same  form  as  Den  Hartog* s 
equation  (3).  Thus  when  b^/b^  m  1'  the  ratio  v^/p^  is  independent 
of  c  and  has  the  invariant  points  property  described  above  in 
conjunction  with  the  conventional  dynamic  vibration  absorber.  The 
value  of  c  that  gives  a  zero  slope  of  the  amplitude  I  v^/F^  *  with 
respect  to  frequency  U)  at  the  invariant  value  of  tVj/F^  I  is  an 
optimum  value  for  c.  Henney  and  Raney  [12]  applied  a  similar 
technique  to  optimize  damping  for  vibrating  uniform  beams  and 
studied  the  sensitivity  of  maximum  displacement  response  to 
deviation  from  optimum  damping. 


Snowdon  [13]  ,[14]  described  how  absorbers  could  reduce  the 
force  transmitted  to  the  terminations  of  undamped  cantilever  beams  at 
the  resonant  frequencies  of  beam  vibration.  The  systems  under 
consideration  are  shown  in  Figure  7.  The  force  transmissibility  T 
across  the  beam  is  defined  as  follows 


T 


»F  /F  I 
T  o 


(20) 


where  F^  is  the  force  transmitted  to  the  fixed  end  and  is  the  applied 
force.  When  a  single  absorber  is  used,  as  in  Figures  7a  and  7b. 
essentially  the  same  procedure  as  that  used  for  the  conventional 
absorber  is  applicable.  The  fixed  points  are  located  by  comparing 
transmissibility  curves  obtained  for  zero  and  infinite  damping. 

The  frequency  ratio  is  then  chosen  so  that  the  fixed  points  lie  on 
transmissibility  curve  actually  take  the  equal  values  of 
transmissibility  at  the  fixed  points.  When  two  absorbers  are  applied 
simultaneously  to  the  beam.  Figure  7c,  Snowdon's  approach  is  to 
assign  each  absorber  the  values  of  optimum  tuning  and  damping  that 
were  determined  when  the  absorber  was  attached,  individually,  to  its 
present  position  on  the  beam.  Optimum  supression  of  both  the  first 
and  the  second  beam  resonances  results  from  this  procedure. 

Jacquot  and  Foster  [15]  considered  an  undamped  single-degree- 
of-freedom  system  equipped  with  a  damped  cantilever  beam  serving  as 
a  vibration  absorber.  The  authors  developed  an  approximate  system 
dynamics  by  an  assumed  modes  approach  retaining  only  one  mode  of  the 
beam.  Since  the  resulting  equations  were  of  the  form  obtained  for 
the  conventional  absorber  problem,  the  same  tuning  and  damping 
approach  was  applicable.  Snowdon  [16]  analyzed  platelike  absorbers 
attached  either  to  a  lumped  mass-spring  system  or  to  a  plate  with 
small  internal  damping  and  presented  optimum  design  parameters  in 
graphical  form. 

Jacquot  [17]  appended  a  discrete  absorber  with  structural  or 
viscous  damping  to  a  beam.  Figure  8.  Using  an  assumed  mode  approach 
as  in  [ 15 ) ,  he  derived  an  approximate  frequency  response  given  by 


(23) 


Because  the  form  of  equation  (21)  is  the  same  as  that  found  in  the 
conventional  absorber  problem  with  the  exception  that  /i  is  replaced 
by  /ip  ( a},  Jacauot  was  able  to  tune  and  damp  to  flatten  the 
fundamental  resonance  response  using  Den  Hartog's  procedure. 

Warburton  and  Ayorinde  [18]  extended  Jacauot' s  representative 
mode  method  to  plates  and  shells  by  defining  appropriate  effective 
mass  and  stiffness  for  such  elastic  bodies.  The  accuracy  of  this 
single-mode  approximation  was  favorably  affected  if  adjacent 
resonant  frequencies  were  well  separated  from  the  natural  frequency 
for  which  the  absorber  was  being  tuned  and  adversely  affected  as  the 
absorber  size  increased.  In  a  companion  paper  [19],  the  authors 
considered  cylindrical  shells  as  examples  of  dynamically  complex 
structures,  for  which  the  ratio  of  adjacent  natural  frequencies 
tends  toward  unity.  They  showed  that  as  dynamic  complexity 
increased  optimum  absorber  parameters  deviated  from  those  calculated 
for  an  equivalent  sing le-degree-of- freedom  system. 

MULT I PARAMETER  DESIGN  BY  NONLINEAR  PROGRAMMING 


McMunn  and  Plunkett  [20]  developed  a  computational  method  to 
optimize  multiple  dampers  for  large  mechanical  systems.  Damping  was 
defined  to  be  optimum  if  the  maximum  response  over  a  range  of 
excitation  frequencies  was  minimized.  The  response  function  of 
interest ,  for  example 


f  (U>,c) 


(24) 


where  x .  is  the  displacement  of  the  ith  mass  and  is  the  force  on 
the  jtn  mass,  is  first  maximized  with  respect  to  excitation 
frequency  u  using 
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and  then  minimized  with  respect  to  the  vector  c  of  damping  values 


min  f (w,c)  (26) 

£ 

In  fact,  the  minimization  over  c  also  considers  values  of  f 
at  u  -  0  ,  and  w  •  fl  if  a  finite  frequency  interval  [0,0]  is  ^^9 
considered.  Two  multlple-degree-of-freedom,  multiple-damper 
discrete  systems  and  a  column  with  complex  modulus  damping  were 
studied  as  examples  using  this  approach.  A  similar  approach  was 
adopted  by  Ng  and  Cunniff  [21]  who  designed  a  three 
degree-of-freedom  isolation  system  and  verified  their  results  by 


experimental  tests.  In  their  formulation,  the  authors  define  a 
primitive  function 

I*!2  -  (P_1f ) (P_1f )*  (27) 

where  *  denotes  complex  conjugate  transpose, 

2 

P  *  -w  M  +  juiC  +  K  (28) 

and  f  is  the  force  vector.  Their  first  step  is  to  minimize 

over  damping  variables  the  primitive  function  maximized  over  frequency 

2 

min  max  I x ■ I  ( 29 ) 

c  CJT 

They  repeat  this  procedure  until  an  optimum  is  reached 

1*1  2  .  *  min  max  I  x°  I  2  (30) 

oot  n  n 
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Dale  and  Cohen  [22]  extended  the  method  of  McMunn  and  plunkett 
to  continuous  systems  whose  steady-state  equations  of  motion  could  be 
reduced  to  a  set  of  ordinary  differential  equations  containing 
spatial  coordinates  as  independent  variables.  Both  dissipative  and 
nondissipative  design  parameters  were  included. 

Lunden  [23]  presented  a  nonlinear  programming  solution  based  on 
a  sequential  unconstrained  minimization  technique  to  the  problem  of 
determining  a  continuous  damping  distribution  which  minimizes  the 
maximum  response  of  a  vibrating  beam  over  a  specified  frequency 
interval.  In  a  second  paper  [24],  the  author  applied  the  same 
approach  to  vibrating  frames.  In  these  references,  the  maximum 
response  F  in  the  frequency  interval  studied  is  written  as 

F  (7?d,7jS)  -  f(T)d,if8,w)  (31) 

d  s 

where  7}  denotes  distributed  structural  damping  and  y  denotes  the 

structural  damping  constants  for  discrete  springs  in  the 

system.  An  exact  displacement  method  is  used  with  hysteretic 

damping  introduced  by  the  loss  factor  7)6  giving  a  complex  bending 

stiffness  EI(  1  +  ii)4).  The  4x4  stiffness  matrix  for  a  beam 

element  then  takes  the  form 

E  -  EI(1  +  x  (32) 

where  K  is  a  4  x  4  matrix  of  transcendental  functions  of  frequency. 


JCitis  [25]  utilizes  structural  reanalysis  and  modal  techniques 
with  nonlinear  programming  to  make  tractable  problems  in  which  the 
systems  under  consideration  contain  a  large  number  of  degrees  of 
freedom.  The  repetitive  computations  of  response  required  in  the 
nonlinear  programming  portion  of  the  optimal design  are  carried  out 
using  efficient  reanalysis  methods  or  condensed  eigenproblem 
solutions  so  that  computation  time  in  the  structural  analysis  phase 
of  the  design  is  reduced. 

IMPEDANCE  MATCHING 

The  progressive  wave  solution  to  the  wave  equation  for 
continuous  chain-like  systems  has  been  used  to  minimize  vibratory 
response  over  a  frequency  range  by  means  of  "impedance  matching". 

To  illustrate  the  basic  idea  of  this  technique  consider  the  rod 
shown  in  Figure  9,  [141.  Designating  the  displacement  at  any  point 
x  by  {,  the  wave  equation  may  be  written  as 
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If  the  end  of  the  rod  at  x  «  o  is  subjected  to  a  sinusoidally 
varying  force  F^,  then  the  driving  point  impendance  Z^  is  given  by 
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Snowdon  [14]  shows  that  the  impedance  Z  defined  in  equation 
(351,  may  be  written  as 


z0  =  aJ?e 


/  - 

I  v-  ReJ<f> 


where  Re^  describes  the  relative  magnitude  of  and  the  phase 
difference  between  the  incident  and  the  reflected  waves.  The 
characteristic  impedance  Z  .  is  defined  as  the  impedance  of  an 
infinitely  long  rod  in  which  reflections  do  not  occur,  the  value  of 

2  w  is  obtained  from  equation  ( 36 )  by  equating  R  to  zero 
ch 
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A  matched  condition  will  occur  if  a  damper  of  damping  constant  Z  ^ 
is  attached  to  the  rod  at  x  ■  0,  This  attachment  causes  the  ratio  R 
of  reflected  to  incident  wave  to  be  zero  so  that  no  vibration 
response  buildup  due  to  reflected  waves  is  possible.  This  idea  has 
been  applied  to  shafts  on  supports  [26].  Here  the  dynamic  response  is 
expressed  in  progressive  wave  form  like  electrical  response  waves  in 
transmisison  line  theory.  Waveforms  for  a  uniform  shaft  flexibly 
supported  on  two  rotational  and  translational  mass-spring-damper 
units  at  the  ends  and  one  such  unit  in  the  interior  are  obtained. 

The  terminating  impedance  is  made  equal  to  the  characteristic 
impedance  of  the  shaft  to  obtain  a  matched  condition. 

LIMITING  PERFORMANCE  METHODS 

A  limiting  performance  approach  was  applied  in  Refs.  [27]  and 
[28]  to  the  optimal  design  of  vibratory  systems  over  a  frequency 
range.  In  this  approach,  instead  of  fixing  the  design  configuration 
at  the  outset,  those  parts  of  the  system  to  be  designed  are  replaced 
by  control  forces.  Then,  for  a  selected  cost  function  and  design 
constraint  the  absolute  optimal  performance  of  the  system  is 
computed  by  solving  an  optimization  problem  in  which  the  control 
forces  axe  unknowns.  The  solution  is  called  the  limiting 
performance  of  the  system.  After  the  limiting  performance 
characteristics  have  been  found,  the  designer  can  choose  a 
prospective  configuration  for  the  part  of  the  system  to  be  designed 
and  apply  parameter  identification  techniques  to  find  optimum  design 
variable  values  so  that  the  designed  system  responds  as  closely  as 
possible  to  the  limiting  performance  response.  The  limiting 
performance  characteristics  are  found  by  linear  programming  and 
parameter  identification  can  be  accomplished  by  such  curve  fitting 
techniques  as  least  squares.  This  two-stage  procedure  has  been 
demonstrated  for  the  optimal  design  of  rotor  suspension  systems  in 
[27]. 

BOOKS  AND  MONOGRAPHS 

The  book  by  Haug  and  Arora  [10]  contains  considerable  material 
applicable  to  frequency  response  shaping,  although  specialized 
aspects  are  not  treated  in  detail.  Two  other  useful  references  are 
the  monograph  by  Sevin  and  Pilkey  [28]  and  the  book  by  Snowdon  [14). 
Sevin  and  Pilkey  present  an  introduction  to  the  subject  and  a 
summary  of  the  state-of-the  art  up  to  1971.  Snowdon’s  book  contains 
a  wealth  of  information  on  vibration  absorbers  and  reduction  of  beam 
vibrations . 
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OJ  '  Wo  =  M  +M  9  Kl  =  spring  constant  of  main 
o  12 

system,  M1  =  mass  of  main  system,  M2  =  mass  of 
absorber,  w  =  excitation  frequency. 

Figure  3.  Transmissibility  of  the  Three-Element 
Dynamic  Absorber 


Figure  4. 


Transmissibility  of  the  Dual 
Absorber 


Figure  5.  Vibration  Absorber  for 
a  Damped  System 


c.  Attached  at  the  end  and  midpoint 


K  =  absorber  spring  constant 
a 

r\  -  absorber  damping  constant 

3l 

M  38  absorber  mass 
a 

Z  «  absorber  impedance 

Figure  7.  Dynamic  Absorbers  Attached  to  a  Cantilever 

Beam  Excited  by  a  Sinusoidally  Varying  Force 
at  its  Free  End 


Figure  8.  Euler-Bernoulli  Beam  with  Dynamic 
Vibration  Absorber 


Figure  9.  Rod  Under  Longitudinal 
Vibration 


